The Annals of Applied Probability 
2009, Vol. 19, No. 6, 2270-2300 
DOI: 10.1214/09-AAP610 

© Institute of Mathematical Statistics, 2009 

LIMIT THEOREMS FOR ADDITIVE FUNCTIONALS 
OF A MARKOV CHAIN 

By Milton Jara\ Tomasz Komorowski^'^ and Stefano Olla^'^ 

CEREMADE Universite de Paris Dauphine, Institute of Mathematics, 
UMCS and CEREMADE Universite de Paris Dauphine 

Consider a Markov chain {X„}„>o with an ergodic probabihty 
measure tt. Let ^ be a function on the state space of the chain, 
with a-tails with respect to tt, a € (0,2). We find sufficient condi- 
tions on the probabihty transition to prove convergence in law of 
N^^°''Y^'^ '^{Xn) to an a-stable law. A "martingale approximation" 
approach and a "coupling" approach give two different sets of condi- 
tions. We extend these results to continuous time Markov jump pro- 
cesses Xt, whose skeleton chain satisfies our assumptions. If waiting 
times between jumps have finite expectation, we prove convergence of 
AT-i/" f^^ V{X,) ds to a stable process. The result is applied to show 
that an appropriately scaled limit of solutions of a linear Boltzman 
equation is a solution of the fractional diffusion equation. 

1. Introduction. Super diffusive transport of energy is generically ob- 
served in a certain class of one-dimensional systems. This can be seen nu- 
merically in chains of anharmonic oscillators of the Fermi-Pasta-Ulam type 
and experimentally in carbon nanotubes (see [20] for a physical review). 
The nature of the stochastic processes describing these emerging macro- 
scopic behaviors is a subject of a vivid debate in the physical literature and 
remarkably few mathematical results are present for deterministic micro- 
scopic models. 
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The macroscopic behavior of the energy in a chain of harmonic oscillators 
with the Hamiltonian dynamics perturbed by stochastic terms conserving 
energy and momentum has been studied in [2]. The density of energy dis- 
tribution over spatial and momentum variables, obtained there in a proper 
kinetic limit, satisfies a linear phonon Boltzmann equation, 

(1.1) dtu{t,x,k)+uj'{k)dxu(t,x,k) = J R{k,k'){u{t,x,k') — u(t,x,k))dk'. 

As we have already mentioned, u{t, x, k) is the density at time t of energy 
of waves of Fourier's mode k € [0, 1], and the velocity uj'{k) is the derivative 
of the dispersion relation of the lattice. 

We remark at this point that (1.1) appears also as a limit of scaled wave, 
or Schrodinger equations in a random medium with fast oscillating coeffi- 
cients and initial data. It is sometimes called, in that context, the radiative 
transport equation (see, e.g., [1, 10, 12, 21, 27], or monography [13] for more 
details on this subject). 

Since the kernel R{k,k') appearing in (1.1) is positive, this equation has 
an easy probabilistic interpretation as a forward equation for the evolution of 
the density of a Markov process {Y{t),K(t)) on M x [0, 1]. In fact, here K(t) 
is an autonomous jump process on [0, 1] with jump rate R{k, k'), and Y{t) = 
Jq uj' {K (s)) ds is an additive functional of K{t). Momentum conservation in 
the microscopic model imposes a very slow jump rate for small k: R{k, k') ~ 
k"^ as /c ~ 0, while velocity uj'{k) remains of order 1 even for small k. So 
when Kit) has a small value, it may stay unchanged for a long time, as does 
the velocity of Y{t). This is the mechanism that generates on a macroscopic 
scale the super diffusive behavior of Y[t). 

The above example has motivated us to study the following general ques- 
tion. Consider a Markov chain {X„,n > 0} taking values in a general Pol- 
ish metric space {E^d). Suppose that vr is a stationary and ergodic prob- 
ability Borel measure for this chain. Consider a function ^ -.E and 
Sn '■= J2n=o ^i^n)- If ^ is centered with respect to vr, and possesses a sec- 
ond moment, one expects that the central limit theorem holds for N~^/'^Sn, 
as N ^ +00. This, of course, requires some assumptions about the rate of 
the decay of correlations of the chain, as well as hypotheses about its dy- 
namics. If ^ has an infinite second moment and its tails satisfy a power 
law, then one expects, again under some assumption on the transition prob- 
abilities, convergence of the laws of N~^/°'S]\f, for an appropriate a to the 
corresponding stable law. 

In 1937 W. Doeblin himself looked at this natural question in his seminal 
article [7]. In the final lines of this paper, he observes that the method 
of dividing the sum into independent blocks, used in the paper to show 
the central limit theorem for countable Markov chains, can be used also 
in the infinite variance situation. A more complete proof, along the line of 
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Doeblin's idea, can be found in an early paper of Nagaev [24], assuming a 
strong Doehlin condition. 

Starting from the early sixties, another, more analytical approach, has 
been developed for proving central limit theorems for Markov chains, based 
on a martingale approximation of the additive functional. By solving (or by 
approximating the solution of) the Poisson equation (I — P)it = ^ where P 
is the transition probability matrix, one can decompose the sum Sn into a 
martingale plus a negligible term, thus reducing the problem to a central 
limit theorem for martingales. This is exploited by Gordin (see [15]) when 
P has a spectral gap. In the following decades, much progress has been 
achieved using this approach. It has found applications in stochastic homog- 
enization, random walks in random environments and interacting particle 
systems (i.e., infinite-dimensional problems, where renewal arguments can- 
not be applied), culminating in the seminal paper of Kipnis and Varadhan 
[18] where reversibility of the chain is exploited in an optimal way (see also 
[5, 6, 14]). For nonreversible chains there are still open problems (see [22] 
and the review paper [25] for a more detailed list). 

As far as we know, the martingale approximation approach has not been 
developed in the case of convergence to stable laws of functionals of Markov 
chains, even though corresponding theorems of martingale convergence have 
been available for a while (cf., e.g., [3, 9]). The present article is a first step 
in this direction. 

More precisely, we are concerned with the limiting behavior of functionals 
formed over functions ^ with heavy tails that satisfy a power law decay, that 
is, 7r(^' > A) ~c+A-" and 7r(^' < -A) ~c-A-" for A » 1 with a G (0,2). We 
prove sufficient conditions under which the laws of the functionals of the 
form N-^/'^Sn converge weakly to a-stable laws, as — > +oo. Theorem 2.4 
is proven by martingale approximation, under a spectral gap condition. 

We also give a proof by a more classical renewal method based on a cou- 
pling technique inspired by [4] . The coupling argument gives a simpler proof 
but under more restrictive assumptions on the form of the probability tran- 
sition (cf. Condition 2.5). We point out, however, that such hypotheses are 
of local nature, in the sense that they involve only the behavior of the pro- 
cess around the singularity. In particular, the spectral gap condition (which 
is a global condition) can be relaxed in this coupling approach, to a moment 
bound for some regeneration times associated to the process (cf. Theorem 
2.7). 

Next, we apply these results to a continuous time Markov jump process 
{A(,t > 0} whose skeleton chain satisfies the assumptions made in the re- 
spective parts of Theorem 2.4. We prove that if the mean waiting time t{x) 
has a finite moment with respect to the invariant measure vr and the tails 
of V{x)t{x) obey the power laws, as above, then finite-dimensional distribu- 
tions of the scaled functional of the form N^^/°' Jq V{Xs)ds converge to the 
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respective finite-dimensional distribution of a stable process (see Theorem 
2.8). 

Finally, these results are applied to deal with the limiting behavior of 
the solution u{t,x,k) of the linear Boltzmann equation (1.1) in the spatial 
dimension d=l. We prove that the long-time, large-scale limit of solutions 
of such an equation converges to the solution of the fractional heat equation 

corresponding to a stable process with exponent a = 3/2. Both approaches 
(i.e., martingale approximation and coupling) apply to this example. 

Note added to the second version: After completing the first version of 
the present paper [16], we have received a preprint by Mellet, Mischler and 
Mouhot [23] that contains a completely analytical proof of the convergence 
of the solution of a linear Boltzmann equation to a fractional diffusion. The 
conditions assumed in [23] imply the same spectral gap condition as in our 
Theorem 2.4; consequently the corresponding result in [23] is related to our 
Theorem 2.8. 



2. Preliminaries and statements of the main results. 



2.1. Some preliminaries on stable laws. In this paper we shall consider 
three types of stable laws. When a G (0,1), we say that X is distributed 
according to a stable law of type I if its characteristic function is of the form 
jggt^x _ gi^CO ^ where the Levy exponent equals 

(2.1) H0-=<^ [ (e^^«-l)|A|-i-"c,(A)dA 

and 

(2 2) c (A) -=1^*' "^^^^ ^ ^ ^' 

* \ct , when A > 0, 

where c~, c+ > and c~ + c+ > 0. The stable law is of type II if a G (1,2) 
and its Levy exponent equals 



(2.3) HO-=a / (e*^«-l-iAe)|A|-i~"c,(A)(iA. 

JR 

Finally, the stable law is of type III is q = 1 and its Levy exponent equals 

(2.4) m-= I (e*«^-l-ieAl[„i,i](A))|A|-2c,(A)dA. 



We say that {Z{t),t > 0} is a stable process of type I (resp., II, or III) if 
Z(0) = and it is a process with independent increments such that Z{\) is 
distributed according to a stable law of type I (resp., II, or III). 
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2.2. A Markov chain. Let {E,d) be a Polish metric space, £ its Borel 
cr-algebra. Assume that {Xn,n > 0} is a Markov chain with the state space 
E and vr — the law of Xq — is an invariant and ergodic measure for the chain. 
Denote by P the transition operator corresponding to the chain. Since vr is 
invariant it can be defined, as a positivity preserving linear contraction, on 
any L^(vr) space for p € [1, +oo]. 

Condition 2.1. Suppose that ^' : ^ M is Borel measurable such that 
there exist a G (0,2) and two constants satisfying c;^ + c~ > and 

lim A"7r(^> A) =c+, 

A— »+oo 

(2.5) 

lim A"7r(^'<-A) = c;. 

A— >+oo 

Condition (2.5) guarantees that ^' E -^^^(vr) for any fi < a. 

In the case of a € (1,2), we will always assume that / \I'd7r = 0. We are 
interested in the asymptotic behavior of Sn := J2n=i ^i^n)- We are look- 
ing for sufficient conditions on the chain, which guarantee that the laws of 
converge to a a-stable law, as ^ +oo. 

We present two different approaches (by martingale approximation and 
by coupling) with two separate set of conditions. 

2.3. The martingale approach result. We suppose that the chain satisfies: 

Condition 2.2 (Spectral Gap Condition). 

(2.6) sup[||P/||i2(,) : / ± 1, = 1] = a < 1. 

Since P is also a contraction in ^-'^(Tr) and L°°(tt) we conclude, via the 
Riesz-Thorin interpolation theorem, that for any p € [1, +oo), 

(2.7) ll^'/llL.W<ai~l'/^"'l||/llL.W, 

for all / G LP{tt), such that / /dyr = 0. 

In addition, we assume that the tails of ^ under the invariant measure do 
not differ very much from those with respect to the transition probabilities. 
Namely, we suppose: 

Condition 2.3. There exists a measurable family of Borel measures 
Q{x,dy) and a measurable, nonnegative function p{x,y) such that 

(2.8) P{x,dy) = p{x,y)Tr{dy) + Q{x,dy) for all x € E', 

(2.9) C{2) :=sup [ p'^{x,y)TT{dx) <+oo 

yeEJ 
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and 

(2.10) Q(x,|^'| > A) <C / p{x,y)TT{dy) yx&E,X>0. 

Jl\'f{y)\>M 

A simple consequence of (2.8) and the fact that vr is invariant is that 

(2.11) J p{x,y)TT{dy) <1 and J p{y,x)TT{dy) <1 Mx^E. 

If a G (1,2) then, in particular, ^ possesses the first absolute moment. 

Theorem 2.4. We assume here Conditions 2.1-2.3. 

(i) Suppose a € (1,2), ^ is centered. Furthermore, assume that for some 
a' > a, we have 

(2.12) ||P^||^.,(^)<+oo. 

Then the law of N-^/'^Sn converges weakly, as N ^ +oo, to a stable law of 
type II. 

(ii) If a (0,1), then the law of N^^^^Sn converges weakly, as N ^ 
+00, to a stable law of type I. 

(iii) When a = l, assume that for some a' > 1, we have 

(2.13) SUp||P^'Ar||^c'f^) < 

N>1 ^ ' 

where ■= < N]. Let cn ■= J^Ndn. Then, the law of N^'^{Sn - 

Ncm) converges weakly, as N ^ +oo, to a stable law of type III. 

Remark. A simple calculation shows that in case (iii) cn = (c+o(l)) log N 
for some constant c. 

2.4. The coupling approach results. 

Condition 2.5. There exists a measurable function 6:E ^ [0,1], a 
probability q and a transition probability Qi{x,dy), such that 

P{x,dy) = e{x)q{dy) + (1 - e{x))Qi{x,dy). 

Furthermore, we assume that 

(2.14) 6:= j e{x)TT{dx)>Q 

and that the tails of distribution of ^ with respect to Qi{x, dy) are uniformly 
lighter than its tails with respect to q, 

t^,.. y Ql(^,|^l>A) „ 

2.15) hm sup 77— -T — = 0. 

^ ' A^oo^gS g(|^|>A) 
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Clearly, because of (2.15), the function ^ satisfies condition (2.5) also 
with respect to the measure q, but with different constants. 

lim A"g(^' > A) = c+^-\ 

A^+oo 

(2.16) 

lim X'^qi^ <-X)=c-e~\ 

A— >+oo 

The purpose of Condition 2.5 is that it permits to define a Markov chain 
{{Xn,dn),n > 0} on £^ X {0, 1} such that 



(2.17) 







= 0\Xn 


= X, 


5n = e) 


= e{x), 








= l\Xn 


= x, 


5n = e) 


= i-e\ 






G A\5n+1 


= 0,Xn 


= x, 


5n = e) 






nxn+1 


e A\5n+1 


= 1, Xn 


= x, 


5n = e) 


= Qi{x, 





for n > 0. We call this Markov chain the basic coupling. It is clear that the 
marginal chain {X„,n > 0} has probability transition P. The dynamics of 
{{Xn,Sn),n > 0} are easy to understand. When X„ = x, we choose Xn+i 
according to the distribution q{dy) with probability 6{x), and according to 
the distribution Qi{x,dy) with probability 1 — 6{x). 

Let Kn be the nth zero in the sequence {5n, n>0}. In a more precise way, 
define kq '■= 0, and for i > 1, 

Ki :=inf{n > Ki_i,(5„ = 0}. 

Notice that the sequence {kj+i — Ki,i> 1} is i.i.d., and E(Kj+i — Kj) = 9~^. 
We call the sequence {Knjn > 1} the regeneration times. 

Observe that, for any i > 1, the distribution of X^. is given by q{dy). In 
particular, is independent of {Xq, . . . , X^j-i}. Therefore, the blocks 

{{Xm^^Ki), • • • , (-'^Ki + l-lj'^Ki + l-l)} 

are independent. The dynamics for each one of these blocks is easy to un- 
derstand. Start a Markov chain {X\,n > 0} with initial distribution q{dy) 
and transition probability Q\{x,dy). At each step n, we stop the chain with 
probability 6{X}^). We call the corresponding stopping time Ati. Each one of 
the blocks, except for the first one, has a distribution {(Xq , 0), [X\, 1), • • • , 
{X^^_^, 1)}. The first block is constructed in the same way, but starts from 
Xq = Xq instead of with the law q{dy). Now we are ready to state: 

Condition 2.6. 

oo 

^n^+"supP(Ki >n|Xo = x) < +oo. 
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Theorem 2.7. Suppose that a € (1,2) and ^ is centered under ir, or 
a € (0,1). Then under Conditions 2.1, 2.5 and 2.6, the law of N^^/'^Sn 
converges to an a-stable law. 

2.5. An additive functional of a continuous time jump process. Suppose 
that {Tn, n > 0} is a sequence of i.i.d. random variables, independent of .7^ := 
cr(Xo,Xi, . . .) and such that tq has exponential distribution with parameter 
1. Suppose that t:E ^ (0,+oo) is a measurable function such that t{x) > 
t^. > 0, X e E. Let 

N 

(2.18) := ^ t(X„)r„. 

n=0 

One can define a compound Poisson process Xt = X^, t G [tiy,tN+i)- It is 
Markovian; see, for example. Section 2 of Appendix 1, pages 314-321, of [17] 
with the generator 

(2.19) Lf{x) = t-\x) j[f{y)-f{x)\P{x,dy), feBb{E). 

Here Bb{E) is the space of bounded and Borel measurable functions on E. 
Let 

(2.20) t:= J tdn <+oo. 

Suppose y — 5- M is measurable and ^'(x) := V{x)t{x) satisfies condition 
(2.5). We shall be concerned with the limit of the scaled processes, 

(2.21) y^(t):=_^y^ V{X{s))ds, t>0, 

as — > +00. Then t~^t{x)TT{dx) is an ergodic, invariant probability measure 
for {Xf,t > 0}. Our result can be formulated as follows. 

Theorem 2.8. (i) Suppose that a G (1,2) and that the assumptions of 
either part (i) Theorem 2.4, or of Theorem 2.7, hold. Then, the convergence 
of finite- dimensional distributions takes place to a stable process of type II. 

(ii) In case a £ (0, 1), we suppose that the assumptions of either part (ii) 
of Theorem 2.4, or of Theorem 2.7 hold. Then the finite distributions of pro- 
cesses {Yfyf{t),t > 0} converge, as N ^-\-oo, to the respective distributions 
of a stable process of type I. 

(iii) When a = 1 and the assumptions of part (iii) of Theorem 2.4 hold, the 
finite distributions of processes {lAr(t) — CNt,t > 0} converge, as N ^ +oo, 
to the respective distributions of a stable process of type III. Here cm := 

J\^\<N^dTT. 
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3. An application: SuperdifFusion of energy in a lattice dynamics. In 

[2] it is proven that the Wigner distribution associated with the energy of 
a system of interacting oscillators with momentum and energy conserving 
noise converges, in an appropriate kinetic limit, to the solution u{t,x,k) of 
the linear kinetic equation 

. ( dtu{t, X, k) + ui'{k) dxu{t, x, k) = Cu{t, x, k), 

\u{0,x,k) = UQ{x,k), 

where {t,x,k) E [0, +cx3) x M*^ x T'^ and the initial condition uq{-,-) is a func- 
tion of class Ci'°(]R'^ X T'^). Here T is the one-dimensional circle, under- 
stood as the interval [—1/2,1/2] with identified endpoints, and T"^ is the 
d-dimensional torus. The function u>{k) is the dispersion relation of the lat- 
tice and it is assumed that Lo{—k) = uj{k) and u{k) ~ \k\ for \k\ ~ {acoustic 
dispersion). The scattering operator C, acting in (3.1) on variable k, is usu- 
ally an integral operator that is a generator of a certain jump process. 
In the case of d= 1, the scattering operator is given by 

(3.2) Cf{k) = I Rik, k')[f{k') - f{k)] dk' 

with the scattering kernel, 

Rik.k') = ^\2s\n^i2TTk)Wi-Kk') 

(3.3) 

+ 2 sin^ (27r/c') sin^ (vr/c) - sin^ (27rA;) sin^ {2-nk')]. 

We shall assume that the dispersion relation is given by a function uj-.T ^ 
[0, +oo), that satisfies a; G C^(T \ {0}) and 

(3.4) Q|sin(7rA;)| < u;(A;) < Cu|sin(7rA;)|, A; G T, 
for some < q < < -l-oo while 

(3.5) lim uj' (k) = :^Crj. 

In the case of a simple one-dimensional lattice, we have uj{k) = c|sin(7rA;)|. 
The total scattering cross section is given by 

(3.6) R{k)= f R{k,k')dk' = -sm^{TTk){l + 2cos'^{'nk)). 

JT 3 

We define t{k) := R{k)~^ since these are the expected waiting times of the 
scattering process. 

Let {Xn,n > 0} be a Markov chain on T whose transition probability 
equals 

P{k,dk') ■.= t{k)R{k,k')dk'. 
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Suppose that {Tn,n > 0} is an i.i.d. sequence of random variables such that 
To is exponentiahy distributed with intensity 1. Let tn '■= t{Xn)Tn, n>0. 
One can represent then the solution of (3.1) with the formula 

(3.7) u{t,x,k)=Euoix{t),k{t)), 

where 

x{t)=x+ / Lo' {k{s)) ds, 



k{t)=Xn, t£[tn,tn+l), 

and k(0) = Xq = k. We shall be concerned in determining the weak limit of 
the finite-dimensional distribution of the scaled process {N-y'^x{Nt),t > 
0}, as — > +00, for an appropriate scaling exponent a. 
It is straightforward to verify that 

(3.8) Tr{dk) = *-^^dk = ^dk, 

R R 

where R := R{k) dk is a stationary and reversible measure for the chain. 
Then P{k,dk') = p{k,k')7r{dk') where 

p{k,k') = Rt{k)R{k,k')t{k') 

and after straightforward calculations, we obtain 

p{k,k') = 6[cos^(7r/c) + cos^(7rA;') — 2cos^(7r/c) cos^(7rA;')] 

X [(1 + 2cos2(7rA;))(l + 2cos'^{Trk'))]~^ 

(3.9) = 6{[|cos(7rA;)| - |cos(7rA;')l]^ 

+ 2|cos(7rA:) cos(7r/c')|[l — |cos(7rA;) cos(7rA;')|]} 

X [(1 + 2cos2(7rA;))(l + 2 cos2(7rA;'))]~^ 

We apply Theorem 2.8 and probabilistic representation (3.7) to describe the 
asymptotic behavior for long times and large spatial scales of solutions of 
the kinetic equation (3.1). The result is contained in the following. 

Theorem 3.1. The finite- dimensional distributions of scaled processes 
{N~'^^'^x{Nt),t > 0} converge weakly to those of a stable process of type II. 
In addition, for any t> 0, x € M, we have 

(3.10) lim [ \u{Nt,N'^/'^x,k)-u{t,x)\'^dk = 0, 

where u{t,x,k) satisfies (3.1) with the initial condition UQ{N~'^/^x,k) , such 
that uq is compactly supported, and u{t, x) is the solution of 

dtu{t,x) = -{-dlfl^u{t,x), 

(3.11) < [ 

u{Q,x) = I uo{x,k) dk. 
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Proof. We start verifying the hypotheses of Theorem 2.8 by finding 
the tails of 

(3.12) ^{k)=uj'{k)t{k) 

under measure vr. Since uj'{k) is both bounded and bounded away from zero, 
the tails of ^{k) under vr are the same as those of t{k). Note that 

(3.13) 7r{k:t{k)>X) = CR\-^/^{l + 0{l)) forA»l, 

and some Cr > 0. This verifies (2.5) with a = 3/2. Since the density of vr 
with respect to the Lebesgue measure is even and ^ is odd, it has a null 
TT- average. □ 

Verification of hypotheses of part (i) of Theorem 2.4- Note that we can 
decompose P{k, dk') as in (2.8) with p{k, k') given by (3.9) and Q{k, dk') = 0. 
Since p{k,k') is bounded, Condition 2.3 and (2.12) are obviously satisfied. 
Operator P is a contraction on L?'{tt), and by the Hilbert-Schmidt theorem 
(see, e.g., Theorem 4, page 247 of [19]) is symmetric and compact. In con- 
sequence, its spectrum is contained in [—1,1] and is discrete, except for a 
possible accumulation point at 0. 

Lemma 3.2. Point 1 is a simple eigenvalue of both P and P^. 



Proof. Suppose 
(3.14) Pf = f. 

We claim that / is either everywhere positive, or everywhere negative. Let 
f~^,f~ be the positive and negative parts of /. Suppose also that /"'' is 
nonzero on a set of positive vr measure. Then f = f^ — f~ and Pf = Pf^ — 
Pf-. Thus /+ = (P/)+ < Yet 

J f^dTi< J Pf+d7r = I f+dTT, 

thus Pf^ = /+. Likewise, Pf~ = f~ . Since for each k we have p{k, k') > 0, 
except for a set of k' of measure vr zero, we conclude that /"*" > vr a.e., 
hence /~ = 0. 

Now we know that PI = 1. We claim that any other / ^ that satisfies 
(3.14) belongs to spanjl}. Otherwise / — cl for some c would suffer change 
of sign. But this contradicts our conclusion reached above so the lemma 
holds for P. The argument for P^ is analogous. 

As a corollary of the above lemma we conclude that condition 2.2 holds. 
Applying part (i) of Theorem 2.8 to iV-2/3 J^^ Lo'{k{s)) ds, we conclude that 
its finite-dimensional distributions converge in law to an a-stable Levy pro- 
cess for a = 3/2. 
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We use the above result to prove (3.10). To abbreviate the notation denote 
Y/v(i) '■= X + N""^/^ uj'{k{s)) ds. Using probabihstic representation for a 
solution of (3.1), we can write 

u{Nt, N^/^x, k) = EkUo{YN{t),k{Nt)) 

(3.15) 

= E / MC,v)^keMi^yN{t) + iijk{Nt)}dC. 

Here uq^^,!]) is the Fourier transform of u{x,k), and is the expectation 
with respect to the path measure corresponding to the momentum process 
{k{t),t > 0} that satisfies fe(0) = k. Since the dynamics of the momentum 
process are reversible with respect to the normalized Lebesgue measure m 
on the torus and is a simple eigenvalue for the generator C, we have 
l!-P*/llL2(m) — > 0, as t — > +00, provided Jj fdk = 0. Suppose that {a]\f,N > 1} 
is an increasing sequence of positive numbers tending to infinity and such 
that aNN~^^'^ — > 0. A simple calculation shows that for any ^,r/ G M and 
eg(x) := e*^^, we have 

\Ek[e^iYNit))er,ikiNt))]-Ek[e^{YNit-taN/N))er,ikiNt))]\^0 

(3.16) 

as ^ +00. 

Using Markov property we can write that the second term under the absolute 
value in the formula above equals 

med^Nit - taN/N))P''^'e^{k{{N - ajv)t))]. 

Let e^(/c) := erj{k) — e^, where := Jjar]ik) dk. By the Cauchy-Schwarz 
inequality, we obtain 

med^Nit - taN/N))P^^'erj{k{{N - aj^)t))] 

(3.17) -EkedYN{t-taN/N))e^\ 

<{Efc|P«^*e,(fc((iV-a^)t))|'}^/l 

The right-hand side of (3.17) tends to in the sense with respect to 
A;GT, as N ^+oo. 

From (3.17) we conclude that 

f f M^,v)medyN{t-taN/N)) 

(3.18) X P''^^e.n{ki{N -aN)t))]dCdk 
-Y. I I uoiC,v)^kedyN{t - taN/N))e^dCdk 

as N ^ +00. Combining this with (3.16), we complete the proof of the 
theorem. □ 
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Verification of hypotheses of Theorem 2. 7. Here we show the convergence 
of N-^l'^x{Nt) by using the couphng approach of Section 4. Define the 
functions 

go (A;) := sin^(27rA;) = 4[sin^(7rA;) - sin'^(7r/c)], 
qi{k) := I sin^(7r/c), 

which are densities with respect to the Lebesgue measure in T. A simple 
computation shows that R{k,k') = 2~'^[qQ{k)qi{k') + qi{k)qo{k')], and there- 
fore, R{k) = 2~'^[(7o(A;) + qi{k)]. The transition probabiHty P{k,dk') can be 
written as 

qo[k)+qi[k) qo[k) + qi{k) 

In particular, in the notation of Section 4, this model satisfies Condition 
2.5 with q{dk') = qo{k') dk' , 9 = qi/{qo + qi) and Qi{k,dk') = qi{k') dk' . No- 
tice that the behavior around of vr and q is the same. Hence, q{'^{k) > 
A) ~ cA~^/^ for A ^ 1. We conclude, therefore, that the function ^'(A:), given 
by (3.12), satisfies (2.16). Observe furthermore that Qi does not depend on 
k and that Qi{k' ,t{k) > A) ~ cA~^/^ for A ^ 1. Due to this last observation, 
condition (2.15) is satisfied. 

We are only left to check Condition 2.6. But this one is also simple, once 
we observe that the sequence {Sn,n > 0} is a Markov chain with transition 
probabilities 

P{dn+1 = l\Sn = 0) = P{6n+1 = 0|5„ = 1) 

"V2 q^(k)qiik)dk 



P{5n+1 = l\5n = 1) 
P{5n+l=Q\5n = ^) 



-1/2 qo{k) + qi{ky 
1/2 ql{k)dk 
-1/2 qQ{k) + qi{ky 
1/2 ql{k)dk 
1/2 qo{k) + qi{k)' 



We conclude that the regeneration time ki satisfies E[exp{7Ki}] < -|-oo for 
7 small enough. Condition 2.6 is therefore a consequence of the fact that the 
transition probability function Qi{k,dk') does not depend on k; therefore, 
we can write 

¥[ki > n\Ko = A:] = (1 - e{k))F[Ki > n - 1]. 
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4. Proof of Theorem 2.7 by coupling. Because of its simplicity, we present 
first the proof of Theorem 2.7 using a basic couphng argument. Let us define 

M{N)=sup{i>0;Ki<N}. 

Note that A4{N) < +oo a.s. An alternative way of defining M{N) is de- 
manding the inequality K;\4(jv) < N < km{N)+i to be satisfied. Then, we 
have 

MiN) 

(4.1) Sn= ^i + ^N, 

1=0 

where 

N 

i=KA1{iV) + l 

In (4.1) we have decomposed Sn into a random sum of i.i.d. random vari- 
ables, {^Pi,i > 1}, and two boundary terms: ipo and Rn. Notice also that 
ktv — f^i is a sum of i.i.d. random variables. Consequently, the law of large 
numbers gives 

(4.2) ^^k = E(k2-^i) and ^^^R-' = 9, 
a.s., as ^ +oo. 

Observe also that when a G (1,2) and ^ is centered, random variable <pi 
is also centered. Indeed, by the ergodic theorem we have that a.s. 

0= lim = lim — - > ifii X — ^ — - =E(j:)i9, 

N^+oc N N^+ooM{N) ^ N ^ ' 

which proves that 

(4.3) Eipi = 0. 

The idea now is that under Conditions 2.5 and 2.6, the random variable 
ipi is equal to ^{X^-) plus a term with lighter tails. Before stating this result, 
we need a simple lemma. 

Lemma 4.1. Let C, he a random variable such that 

lim x"P(C>x) = c+, lim j;"P(C< -x) =c~. 

x^oo x— ►oo 

Let i he such that lim^.^oo IP(|?| > x)/f{\C\ > x) = 0. Then 

(4.4) lim 2;"P(C + C>a;) = c+, lim x"P(C + C < -x) = c". 



x^oo 
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Proof. Without loss of generality, we just consider the first limit, the 
second one follows considering — ^, — ^. We will prove that the lim inf ^^-.^oo 
of the previous expression is bigger than c+ and the limsup is smaller than 
c+. We start with the upper bound: for any e > there exists xq so that for 
x>xq, we have 

a;"F(C + C>x)< x"P(C > (1 - e)x) + a;"P(^ > ex) 
c+ + e F{\C\>€x) ^ c+ + c_ 



(l-2e)° P(|C|>ex) (e/2)" 
Now take above the upper limit, as x — > +oo, to get 

limsupx°P(C + e > x) < j^^^. 

Since e is arbitrary, we have proved the upper bound. The lower bound is 
very similar: 

¥{C + ^ > x) = ¥{C + ^ > x,^ > -ex) +¥{C + ^> x,^< -ex) 
>P(C>(l + e)x,C>-ex) 

> P(C > (1 + e)x) - P(C > (1 + e)x, C < -ex) 

> P(C > (1 + e)x) - P(C < -ex). 

Starting from this last expression, the same computations done for the 
upper bound show that 

liminf x"P(C + ^>x)> 

Since e > is arbitrary, the lemma is proved for the first expression in 
(4.4). The second case can be done in the same fashion. □ 

Lemma 4.2. Let ^ satisfy (2.5) with constants c^,c~ together with Con- 
ditions 2.5 and 2.6. Then the law of each (pi satisfies 

hm A"P(v?,>A) = c+ri, 

A— >+oo 

(4.5) 

lim A"P((/?i < -A) = c:7^^^ 

A— >+oo 

Proof. The idea of the proof is simple. Random variable is the sum 
of a random variable with an a-tail, ^{X^.), and a finite (but random) 
number of random variables with lighter tails ['^{Xf^.j^i), . . . ,'^{X^.^-^-i)). 
By Condition 2.6, the random number can be efficiently controlled. To sim- 
plify the notation, assume that Xq is distributed according to q, so the first 
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block is also distributed like the other ones. Then 

/ K\ — l \ OO /n— 1 



(4.6) 



j=l / n=l \j=l J 

oo n— 1 

n=l j=l 



The probability under the sum appearing in the last expression can be esti- 
mated by 

¥{^{X,)>t/{n-l),Sp = l,yp<j) 

(4.7) 

= E[Qi(X,_i, ^ > t/{n - l)),6p = l,yp<j]. 

When j >n/2 we can use (2.15) to bound the expression on the right-hand 
side of (4.7) from above by 

(4.8) 7/ = 1, Vp < j] < > K > n/2]. 

Here g{x) is a bounded function that goes to 0, as x — > 00. On the other 
hand, when j < n/2, we rewrite the probability appearing under the sum on 
the right-hand side of (4.6) using the Markov property. It equals 

K[K[ki =n-j\Xo = Xj],^{Xj) > t/{n - l),6o = ■ ■ ■ = 6^ = I]. 

Here {Xn,n > 0} is another copy of the Markov chain {Xn,n > 0}, and ki 
is the respective stopping time defined in correspondence to ki. We can 
estimate this expression by 

F[^{Xj) > t/{n - 1)] supP[Ki >n- j\Xo = x] 

X 

<^^^^^supPK>n/2|Xo=rr]. 

t X 

Summarizing, we have shown that the utmost left-hand side of (4.6) can be 
estimated by 

E ^^^^ supP,(Ki > n/2). 

1 I X 

n=l 

We conclude that this expression is o(t~") by invoking Lebesgue dominated 
convergence theorem and Condition 2.6. The negative tails are treated in 
the same way. Therefore, ipo — ^{Xq) has lighter tails than ^(Xq) itself. By 
Lemma 4.1, the sum of a random variable satisfying condition (2.16) and a 
random variable with lighter tails also satisfies condition (2.16) for the same 
constants c'^6~^, c^9~^. 
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At this point we are only left to recall the classical limit theorem for i.i.d. 
random variables. It follows that there exist 

^ 1=0 J=KA1{Ar) + l 

^(JV)V/" 1 , 1 ^ 



^ ^ 1=0 J=K^(jv)+l 

Recall (4.2), and notice that by (4.5), 

lim A"P(^^/"(/?, > A) = c+ 
lim A"P(^^/"(^i < -A) = c-. 

Let C* := (c~ + c+)^^^. By virtue of the stable limit theorem for i.i.d. ran- 
dom variables (see, e.g., [8], Theorem 7.7, page 153), we know that for cat := 
NE[^i, \^i\< (C*Af)^/"] such that the laws of iV"^/"(EiIo O^^'^fi-CN) con- 
verge to an a-stable law. When a <1 constants cat ~ cN^/°' and they can 
be discarded. Observe, however, that since Ey^i = 0, cf. (4.3), for a S (1,2), 
we have 

/• + CXD 

CN = -N¥.[ipi,\^i\> {C^Nf''^]= N / [P[v9i < -A] -F[v9i > A]](iA 



r+oc fJ\ 

-'N (c;c+)- = C7(l + o(l))ivV" 

JfCTV i/« A°^ 



i{c,Ny/° 

for some constant C. The constants ctv can again be discarded. We conclude 
therefore that the laws of 



weakly converge to some a-stable law u^,. Since Cn '■= 0"^^°^N~^Ai{N) con- 
verges a.s. to 1, the joint law of {1Cn-,Cn) converges to 5i, as N ^ +co. 
According to the Skorochod representation theorem there exists a probabil- 
ity space and random variables {ICn,^n) such that {jC]y,C]y) = {K,n-,Cn) 
for each N and {1Cn-,£-n) — > (X*-, 1) a.s. The above in particular implies that 
K-NCjs, converges a.s. to Y^. Since K-nCn ~ ^nCn-> conclude the conver- 
gence of the laws of ICnCn to z/^. □ 
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5. Proof of Theorem 2.4 by martingale approximation. Below we formu- 
late a stable limits law that shall be crucial during the course of the proof 
of the theorems. 

Suppose that {Zn,n > 1} is a stationary sequence that is adapted with re- 
spect to the filtration {Qn, n > 0} and such that for any / bounded and mea- 
surable, the sequence {E[/(Z„)|C/„_i], n > 1} is also stationary. We assume 
furthermore that there exist a G (0, 2) and c;^, c~ > such that + c~ > 
and 



(5.1) 



P[Zi>A]=A-"(c+ + o(l)), 

P[Zi < -A] =A-"(c; + o(l)) as A 
In addition, for any g G Cq°(IR \ {0}), we have 



(5.2) 

and 
(5.3) 



lim E 



[Nt] 
n=l 



at / g{X) 



+ 00. 



c*(A) dX 



|A|i+° 



lim NE<E 



n > 2 



0. 



Here c*(-) appearing in (5.2) is given by (2.2). Let Mn := J2n=i ^n, N >1 
and Mo := 0. 

When a = 1 we shall also consider an array {Zi^^:n>l}, N>1 of sta- 
tionary sequences adapted with respect to the filtration {Qn : n > 0}. Assume 
furthermore that for each A'^ > 1 and any / bounded and measurable se- 
quence, {K[f{Zn^^)\Gn-i] :n > 1} is stationary. We suppose that there exist 
nonnegative c'^ , c~ such that cf + c~ >0 and 

(5.4) lim sup[|AP[zj^^>A]-c+| + |AP[z{^)<-A]-c;|]=0. 

A-»+oo^>l 



Let 



N 



Mm := - IE[^^^1[I^^^I < NWGn^,]}, iV > 1, 



71=1 



and Mq := 0. The following result has been shown in Section 4 of [9]. 

Theorem 5.1. (i) Suppose that a G (1,2), conditions (5.1)-(5.3) hold, 
and 

(5.5) E[Z„|a„,„i]=0 forn>l. 

Then A^~-^/"M[jv.] ^ Z{-), as N ^ +oo, weakly in D[0, +oo), where {Z{t),t > 
0} is an a-stable process of type II. 
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(ii) Suppose that a € (0, 1) and conditions (5.1)-(5.3) hold. Then N^^^" x 
M[7v-] =^ ^i'); as N ^ +00, weakly in D[0,+oo), where {Z{t),t > 0} is an 
a-stable process of type I. 

(iii) For q = 1, assume (5.2) and (5.3) with replacing Z„, and (5.4)- 
Then A^~^M[jv.] =^ as N ^ +oo, weakly in D[0, +oo) to a Levy process 
{Z{t),t>0} of type III. 

Proof of part (i) of Theorem 2.4- Let x ^ L^i'^)^ P ^ (Ij be the unique, 
zero-mean solution of the equation 

(5.6) x-Px = '^- 

Since ^' G L^iT^) for /? S (0,a) is of zero mean, the solution to (5.6) exists 
m L^(7r) and is given by x = J2n>o ^"^^ ■ This follows from the fact that 
< a(^/^~^)"||^'||^/3, n > [see (2.7)], so the series defining x geo- 
metrically converges. Uniqueness is a consequence of (2.7). Indeed, if xi was 
another zero-mean solution to (5.6), then 

IIX - XiWlp = \\P{X - Xi)\\l0 < - XiIIl/3, 

which clearly is possible only when X ~ Xi = (recall that a < 1). Note also 
that from (2.12) it follows that in fact Px={I- P)~^{P^) € L"'(7r). Thus 
in particular, 

(5.7) vr(|Pxl>A)< 

and consequently x satisfies the same tail condition as ^ [cf. (2.5)]. 
Then by using (5.6), we can write 

N N 

(5.8) Sn = Y. ^(^n) = E ^« + ^x(^o) - Px{Xn) 

n=l n=l 
with Zn = x{Xn)-Px{Xn-l). 

In what follows, we denote by Cq°(R \ {0}) the space of all C°° functions 
that are compactly supported in R \ {0}. According to part (i) of Theorem 
5.1, we only need to demonstrate the following. 

Proposition 5.2. For any 5 G C^(M \ {0}), equalities (5.2) and (5.3) 
hold. 



More explicitly, we have 



E 



g{N-'/''[xiy)-PxiXn^i)])PiXn^i,dy) 
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and using the stationarity of vr, we can bound the left-hand side of (5.2) by 

[Nt] 

P{Xn-udy) 



E 



(5.9) 



E J\9{N-'/''[x{y)-Px{Xn^i)])-9(^) 

n=l 



+ 



^c*(A) dX 

so (5.2) is a consequence of the fohowing three lemmas, each taking care of 
the respective term of (5.9): 

Lemma 5.3. 



(5.10) 



lim N U\g{N~^'''[x{y)-Px{x)]) 

iV— >oo 



Lemma 5.4. 

N 

(5.11) e' 



n=l 

Lemma 5.5 



(5.12) lim 



9{N-^/''x{y))\P{x,dy)7Tidx) =0. 



Y,ng{N-'/-x{Xn))\Gn^i]-N J g{N'^'^x{y))<dy) 



0. 



N J g{N-y-x{y)Mdy)-a J^g(xf-li^ 



0, 



where c*(-) is given by (2.2). 



Proof of Lemma 5. 3. Suppose that supp g C [—M, M] \ [—m,m] for some 
<m< M < +00 and 9>0. Denote 

AN,e = {{x,y) : \x{y) - OPx{x)\ > iV^/"m}. 
The left-hand side of (5.10) can be bounded from above by 

iV^-i/- f'de [ [\g\N~y"[xiy) -ePxix)])Pxix)\P{x,dy)7r{dx) 



< CAf^~^/° / dO 



'deff \Px{x)\P{x,dy)7r{dx) 

J J ^N,0 



< CAfi-i/° / de 



l~l/a' 

Pix,dy)7r{dx)] \\Px\ 
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From the tail behavior of x of [see (5.7) and the remark below that 
estimates], it is easy to see that for any 6 G (0, 1), 

P{x,dy)T:{dx) <¥[\x{Xi)\ > (miV^/")/2] 

+ F[e\PxiXo)\>{mN'/'')/2] 

= ^(1 + 0(1)) 
as ^ 1. Since a' > a we obtain (5.10). 

Proof of Lemma 5.4- To simplify the notation we assume that suppflf C 
[m,M] for < m < M < +oo. Denote Bn,x = {y : xiv) > N^^"'^}- We can 
rewrite the left-hand side of (5.11) as 

N 



poo ^' 

/ g'{\)Y,GN{Xn-i,\)^ 



(5.13) E 
where 

Gn{x, A) = P(x, Bn^x) - tt{Bn,x). 
Notice that / Giy{y, A)7r((iy) = and 

G%{y,X)7T{dy)= J P\y,BN,x)7:{dy)-7:\BN,> 



N,X) 



(5.14) -^I{Ib ^(^y) 

+ 2 J Q\y,BN,xMdy)-7r\B_ 

To estimate the first term on the utmost right-hand side, we use the Cauchy- 
Schwarz inequality, while for the second one we apply condition (2.10). For 
A > m, we can bound the expression on the right-hand side of (5.14) by 

7;'^iBN,m) [ [ p^{x,y)TT{dx)TT{dy) + CTT'^{BN,m.) 

(5.15) ^ 

< — o(l) asN^oo, 

by virtue of (2.9) and the remark after (5.7). Thus we have shown that 

(5.16) iVsup / G%iy,XHdy)^0 

X>mJ 
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as — > oo. We will show now that (5.16) and the spectral gap together 
imply that 

2 



(5.17) 



sup E 



N 



GAr(X„_i, A) 



n=l 







as N ^ oo. Since supp^i' C [m, M] expression in (5.13) can be then estimated 
by 



sup E 

X>m 



N 



G'7v(X„_i, A) 



n=l 



X / \g'{\)\d\^^ 



as N ^ +00 and the conclusion of the lemma follows. 

To prove (5.17) let un{-,\) = (/ - -P)~^GAr(-, A). By the spectral gap 
condition (2.7), we have 



(5.18) 



N-l 



l-a2 

We can then rewrite 

N 

GAr(X„_i, A) = U^iXo) - Uiy{XN) + 
n=l n=l 

where Un = u^iXn) — Puj\[{Xn-i), n > 1 is a stationary sequence of mar- 
tingale differences with respect to the natural filtration corresponding to 
{Xn,n>0}. Consequently, 

2 



E 



N 



GN{Xn-l, X) 



n=l 



<CN u%iy,XMdy)^0 



and (5.17) follows from (5.16) and (5.18). 

Proof of Lemma 5.5. To avoid long notation, we again assume that 
supp5( C [m, M] for < m < M < +oo. The proof in the case oig C [— M, —m] 
is virtually the same. Note that 

x{y) 



N 9 



7r{dy) 



N 



7r(x > A) d\ 



N g'{\)7T{x>N'l^\)d\. 
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Thanks to (2.5) the last expression tends, however, as — > +00, to 



c+ dX 



c*(A)dA 



1+1 



Proof of Proposition 5.2. We have already shown (5.2), so only (5.3) 
requnes a proof. To simplify the notation we assume (5 = 0. Suppose that 
supp g C [m, M] for some < m < M. We can write 

Zi 



E 



g{N-^/^^{y))p{Xo,y)T^{dy) 



(5.19) 

+ /i(^o, y)5'(A^~'/"(^(y) + Oh{X^,y))) 

xp{Xo,y)7r{dy)de, 

where h{x,y) := Pxiv) ~ Pxi^)- Denote by Ki and K2 the first and the 
second terms appearing on the right-hand side above. By Cauchy-Schwarz 
inequality, 

E( / |Px(y)b(^o,2/)vr(d?/)'' 



llfl'lP r 



(5.20) 



< 



2\\g'\\UPx\\h 



+ E( / |Px(^o)|p(^o,y)^(c^y) 
(^) 



Hence hmAr^+oo A^EK| = 0. 
On the other hand, 

(5.21) < \\g\\^ J p{Xo,y)l[\^{y)\>mN^/"/2]7r{dy), 

and in consequence, by Jensen's inequality, 

(EK2)1/2 < ii^ii^ f (Ep2(Xo,y))V2i[j^(y)| > aiVV"/2]vr(dy) 



<ll5l 



/(x,y)l[|M/(y)| > aiVi/72]7r(dx)7r(dy) 



1/2 



Thus we have shown that 

NKKf < iV7r[|1'| > aA^^/"/2] 
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(5.22) 

X 11 p\x,y)l[\^iy)\>aN^/''/2]7r{dx)7r{dy)^0. 

Condition (5.3) is then a consequence of (5.20) and (5.22). 

Proof of part (ii) of Theorem 2.4- The proof of this part rehes on part (ii) 
of Theorem 5.1. The following analogue of Proposition 5.2 can be established. 



Proposition 5.6. Suppose that a e (0,1). Then for any g ^ C§ 
{0}), 



(5.23) 

and 
(5.24) 



hm E 



m 

n=l 



lim iVE<!E 



Sn- 



t / 5(A) 



C,iX)dX 



|A| 



l+a 



Go 



Proof. The proof of this proposition is a simplified version of the ar- 
gument used in the proof of Proposition 5.2. The expression in (5.23) can 
be estimated by 



E 



(5.25) 



[Nt] 

E 

n=l 



+ 



m / g 



P{Xn-i,dy)-[Nt] J g 
*(y) 



TT{dy) -at g{X) 



c*(A) dX 



TT{dy) 



|A| 



l+a 



The proof that both the terms of the sum above vanish goes along the 
lines of the proofs of Lemmas 5.4 and 5.5. We can repeat word by word the 
argument used there, replacing this time x by As for the proof of (5.24) it 
is identical with the respective part of the proof of (5.3) (the one concerning 
term Ki). □ 

Proof of part (iii) of Theorem 2.4. Recall that ■= < N]. Let 

Xn be the unique, zero mean solution of the equation 



(5.26) 

We can then write, 



XN - PXN = ^Af - CAT. 



TV N 

(5.27) Sn - NcN = E(*(^") - CTv) = E 4"^^ + ^XTv(^o) - ^'xtv(^tv) 
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with 

Z(^) = XN{Xn) - PXN{Xn~l) + ^ {Xn)l[\^ {Xn)\ > N]. 

We verify first assumptions (5.2), (5.3) and (5.4). 

Condition (5.4) is an obvious consequence of the fact that 

(5.28) = PXNiXn) - PXN{Xn^l) + ^'(^n) " Cn 

and assumption (2.13). To verify the remaining hypotheses, suppose that 
suppg C (m,M) and m < 1 < M. Let us fix (5 > 0, to be further chosen later 
on, such that mKl — 5<l + 5< M . We can then write g = gi + g2 + 93 
where each g^ e C°°(]R), H^iHoo ^ Ibllooi and the supports of gi,g2,g3 are 
correspondingly contained in (m, 1 — 5) , {1 — 5, 1 + 5) , {1 + 5, M). We prove 
(5.2) and (5.3) for each of the function ^jS separately. Note that 



E 



9i 



-in 



giiw^''HXn^^,y))P{Xn^i,dy), 



where 

w^''\x,y) :=N-'^{y)l[\^{y)\>N]+N-'[xN{y)-PXN{x)]. 

For i = 1 and i = 3, we essentially estimate in the same way as in parts (i) 
and (ii) of the proof of the theorem, respectively. We shall only consider here 
the case i = 2. 

Note that then w^^\x,y) = w[^\x,y) where 

(5.29) wf^ {x,y) = N-'[^{y) - cn] + N-'eRN{x, y) 

with RN{x,y) := PxN{y) - Pxn{x). However, 

g2{w'^''\Xn-uy))=g2{N~\^{y)-CN)) 

+ N~^RN{Xn^i,y) f' g'2{wf\Xn^uy)) dd 



and 



E 



N 

E 

n=l 



<E 



52 [w--'-' (X„„i , y))P(X„_i , dy) 
^ j g2{N-\^>{y)-CN))P{Xn-udy) 



n=l' 



+ 



152 



2(w^^^ {x, y))RN{x, y)\P{x, dy)-K{dy) dO. 



Denote the first and the second term on the right-hand side by and 



J2^\ respectively. Term J^" ' can be now estimated as in the proof of part 
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(ii) of the theorem. We conclude then, using the arguments contained in the 
proofs of Lemmas 5.4 and 5.5 that 

"1+5 dX 



hmsup jI^'' < 

N- 



\9\ 



♦ + 00 



l-<5 



On the other hand, to estimate hm7v-»+oo J2^^ = 0, since g2{wg" ' {x , y)) — > 
in measure P{x, dy)iT{dy) dO and the passage to the hmit under he integral 
can be substantiated thanks to (2.13). 

Choosing now sufficiently small 5 > we can argue that the calculation of 
the limit can be reduced to the cases considered for gi and (73 and that con- 
dition (5.2) can be established for The proof of (5.3) can be repeated 
from the argument for part (i) of the theorem. 
Finally, we show that 



A2- 



(5.30) 



lim — E 



N 



^E[Zr)l[|Zr)|<iV]|£?„„i] 



n=l 



0. 



Denote the expression under the limit by L^^\ Let A > 1. We can write 



= LY' + L)^' + L)^' depending on whether |^(X„)| > AiV, |^(X„)| G 
(A-iiV,AiV], or \^{Xn)\ < {Ay^N. Then 



L 



(N) 



N 



< ^ P[|*(X„)| > A7V, < iV] = Arp[|^(Xi)| > AiV, |z| 

n=l 

AN) 



<N]. 



From formula (5.28) for Zn ' , we conclude that the event under the condi- 
tional probability can take place only when \PxN{Xn)\, or \PxNiXn-i) \ > 
A^(A — l)/3 for those N, for which CAr/A^<(A — l)/3. Using this observa- 
tion, (2.13) and Chebyshev's inequality, one can easily see that 



4^)<2iV[iV(A-l)/3]-'||Px^l 



L"'(7r) 







as — > -|-oo. To deal with Lg^^ consider a nonnegative g € C°°(M) such 
that ll^lloo < li [^~"^)^] C supp5( C [Aj""^, Ai] for some Ai > A. Repeating 
the foregoing argument for g2, we conclude that 

'•^i dX 



limsupL2^'* < \\g\\oo 



^1 



1 X^' 



which can be made as small as we wish by choosing Ai sufficiently close to 



1. As for L 



(N) 



(5.31) L 



(TV) 



, note that it equals 

N 



1 



E 



^E[AfWi[|MW| <iv,|^(x„)| < {Ay'N]\g. 



71-1 



n=l 
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where 

M^^"") ■. = XN{Xn) - PXN{Xn-l) 

(5.32) 

= ^N{Xn) -CN + PXN(Xn) " PXN (Xn-l). 

Thanks to the fact that Mn^^ are martingale differences, the expression in 
(5.31) can be written as -^^3^^ = — (Lg^^ + L^l^^ + l'^^) where -^^3^*^^ corre- 
spond to taking the conditional expectation over the events Ai for i = 1, 2, 3 
given by 

^i:=[|AfW|>iV,|^(X„)|<(A)-iiV], 
A2 := [|mW| > N, |^(X„)| > (A)-iiV], 
^3 :=[|MW|<iV,|vf(X„)|>(A)-iiV]. 

To estimate -^^3^^ z = 1,2 we note from (5.32) that |-/Vfn^''| > only when 
^NiXn) = ^'(X„) and |^(X„)| < N, or Px7v(^n-i), PXN{Xn) are greater 
than cN for some c > 0. In the latter two cases we can estimate similarly to 
In the first one, however, we end up with the limit 

1 ^ 

limsup-E^E[(|^^(X„)| + \cn\ + \PXN{Xn)\ + |Px^(^n-i)|), 

N>\^{X^)\>A-^N\gn-i] 
< limsupiV(l + \cn\/N)tt[N > |*| > (A)-iiV] 

+ limsup / {I + P)\PxN\l[N>\'i>\ > {A)-^N]d7r. 

The second term on the utmost right-hand side vanishes thanks to (2.13). 
The first one can be estimated as in the proof of Lemma 5.5, and we obtain 
that it is smaller than C J^~i dX, which can be made as small as we wish 
upon choosing A sufficiently close to 1. We can estimate, therefore, 

lim sup L^^^ 



1 ^ 

<lmisnp-EY,m^N(.Xn)\ + \cN\),N >\^{Xn)\> A-'N\gn-i] 

n=l 



X 



1 ^ 

+ hmsup -E ^ E[(|Px^(^n)| + \PXN{Xn-l)\), 

iV->+oo X 

\^{Xn)\>A-'N\gn^i] 
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:limsupiV7r[Af > 1*1 > (A)~iAf] 



N- 



+ limsup / (I + P)|Px7v|l[|^'| > {A)-^N]d7r 



"1 dX 

/A-1 A2' 

which again can be made arbitrarily small. 



<C 



6. Proof of Theorem 2.8. Suppose that we are given a sequence of i.i.d. 
nonnegative random variables {pn^i^ > 0} independent of > 0} and 

such that Aa ■= J^^"^ p"'(p{dp) < +oo, where (p{-) is the distribuant of po 
and a € (0,2). We consider a slightly more general situation than the one 
presented in Theorem 2.4 by allowing 

[Nt] 

(6.1) SN{t) := '^iXn)Pn- 

n=0 

Observe that, if vr is the law of X„, observable ^ satisfies the tail conditions 
(2.5), and pn is independent of Xn, then 

A"P(^'(Xo)po > A) = / A"7r(* > Xp-'Mdp) ctA^. 

Jo A-++00 

Define also 

(6.2) Cn:= [ ^dirEpQ. 

J\'i\<N 

Consider then the Markov chain {(X„,, /)„), n > 0} on x M^. This Markov 
chain satifies all conditions used in the previous sections, with stationary er- 
godic measure given by 7r{dy) <p{dp). Then with the same arguments as 
used in Section 5 we get the following. 

Theorem 6.1. (i) Under the assumptions of the respective part (i), 

or (ii) of Theorem 2.4, we have N'^^"^ Sn{-) =^ as N ^ +oo where 

{Z{t),t > 0} is an a-stahle process of type either type I, or II with the pa- 
rameters of the corresponding Levy measure [cf. (2.2)] given by 



(6.3) a (A) := 



aAaC^ , when A < 0, 
aAacf, whenX>0. 



Here denotes the convergence in the sense of finite- dimensional distribu- 
tions. 

(ii) In addition, under the assumptions of part (iii) of Theorem 2.4 finite- 
dimensional distributions of N~^SN{t) — Cj\[t converge weakly to those of 
{Z{t),t >0}, a stable process of type III. Here Cn is given by (6.2). 
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Remark. The results of the first part of the above theorem follow under 
the conditions of Theorem 2.7, by using the coupling argument of Section 4. 

Let us consider now the process Yj\i{t) defined by (2.21). We only show 
that one-dimensional distributions of Yi\f{t) converge weakly to the respec- 
tive distribution of a suitable stable process {Z{t),t > 0}. The proof of con- 
vergence of finite-dimensional distributions can be done in the same way. 

Given t > define n{t) as the positive integer, such that 

tn{t) <t< t„(t)+i, 

where is given by (2.18). Let 

s{t) := t/t, 

[Nt] 

i?^(i):=iV-V"^^(Xfc)rfc, t>0, 
k=0 

where, as we recall, ^'(x) := V{x)t{x), x G E and {r^, A; > 0} is a sequence 
of i.i.d. variables distributed according to an exponential distribution with 
parameter 1. Using the ergodic theorem one can easily conclude that 

. . , . n(Nt) 

(6.4) SAf(t) := jy s{t) asA^^+oo, 

a.s. uniformly on intervals of the form [to;^] where < to < ^- We have 

n{Nt)-l Ni--t , . 

k=0 



Note that 



^ n{Nt) 



- J2 ^iXk)rk = B^{sN{t)). 



k=0 



Lemma 6.2. For any t>0 and e > fixed, we have 

(6.5) lim F[\YN{t)-BN{sNm>£]=0. 

A*— >+oo 

Proof. Let cr > be arbitrary. We can write that 

F[\YN{t)-BN{sNm>e] 

(6.6) <n\sN{t)-s{t)\>a] 

+ F[\sN{t)-sit)\<a,\YNit)-BNisNm>^]- 
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The second term on the right-hand side can be estimated from above by 

< P[sup{|^'(Xfc)|rfc : k G [{s{t) - a)N, {s{t) + a)N]} > N^/'^e]. 

Using the stationarity of {\^{Xk)\Tk-,k > 0} the term on the right-hand side 
equals 

P[sup{|^'(Xfe)[rfe : A; G [0, 2(jiV]} > iVi/°e] 

r+oQ (~<fj 

< 2aN J e"^7r[|^(x)| > r-^iV^/^e] dr < — 

for some constant C > 0, by virtue of (2.5). From (6.6) we obtain, therefore, 

hmsupP[|y7v(t) - BNisNim >e]< — 

for an arbitrary a > 0, which in turn implies (6.5). □ 

It suffices, therefore, to prove that the laws of -Bat (sat (t)) converge, as A'" ^ 
-|-oo, to the law of the respective stable process. According to Skorochod's 
embedding theorem, one can find pairs of random elements {BN{-),siy{t)), 
N >1, with values in D[0,+oo) x [0, -|-oo), such that the law of each pair 
is identical with that of {B]\f{-), sj\[{t)), and {B]\f{-), S]y{t)) converges a.s., as 
N — > -|-oo, in the Skorochod topology to {Z{-), s{t)). Here {Z{t),t > 0} is 
the stable process, as in Theorem 6.1. According to Proposition 3.5.3 page 
119 of [11], the above means that for each T > there exists a sequence of 
increasing homeomorphisms Aat: [0,r] — > [0,T] such that 

(6.7) hm 7(Ajv) = 0, 

Ai— >+oo 



where 



and 



j{Xn) ■= sup 

0<s<t<T 



AAr(t) - AAr(s) 

log- 



t 



(6.8) sup \Bn o A7v(i) - Z{t)\ = 0. 

te[o,T] 

As a consequence of (6.7) we have of course that 

(6.9) lim sup |AAr(t) -t| = 0. 

Note that the law of each Bj\f{sN{t)) is identical with that of BN(sN{t)). 
We also have 

\BNisN{t)) - Z{s{t))\ < \BN{sNit)) -ZoX-\sNm 

+ \ZoX~\sN{t))-Z{sm. 
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The right-hand side, however, vanishes N +00, thanks to (6.8), 

(6.9) and the fact that for each fixed s > one has P[Z(s— ) = Z{s)] = 1 
(see, e.g., Theorem 11.1, page 59 of [26]). The above allows us to conclude 
that \B]\f{sN{t)) — Z{s{t))\ a.s., as — > +00, thus the assertions of 
Theorem 2.8 follow. 
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